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Abstract 

Let s and t be vectors of positive integers with the same sum. We study the 
uniform distribution on the space of simple bipartite graphs with degree sequence s 
in one part and t in the other; equivalently, binary matrices with row sums s 
and column sums t. In particular, we find precise formulae for the probabilities 
that a given bipartite graph is edge-disjoint from, a subgraph of, or an induced 
subgraph of a random graph in the class. We also give similar formulae for the 
uniform distribution on the set of simple directed graphs with out-degrees s and in- 
degrees t. In each case, the graphs or digraphs are required to be sufficiently dense, 
with the degrees varying within certain limits, and the subgraphs are required to be 
sufficiently sparse. Previous results were restricted to spaces of sparse graphs. Our 
theorems are based on an enumeration of bipartite graphs avoiding a given set of 
edges, proved by multidimensional complex integration. As a sample application, 
we determine the expected permanent of a random binary matrix with row sums s 
and column sums t. 



This is a revised expanded version of "Asymptotic enumeration of dense 0-1 matrices with specified 
line sums and forbidden positions" , submitted to arXiv on 22 Jan 2007. A concise version will appear in 
Random Structures and Algorithms. 

T Research supported by the Australian Research Council. 
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1 Introduction 



Let s = (s 1; . . . , s m ) and t = (t 1; . . . , t n ) be vectors of positive integers with YljLi s j = 
Sfc=i tk- Define B(s, t) to be the set of simple bipartite graphs with vertices {u\, . . . , u m }U 
{vx, . . . ,v n }, such that vertex Uj has degree s,- for j = 1, . . . , m and vertex v k has degree £ fe 
for k = 1, ... ,11. Equivalently, we may think of B(s, t) as the set of all m x n matrices 
over {0, 1} with jth row sum equal to Sj for j = 1, . . . , m and kth column sum equal to t k 
for k = 1, . . . , n. 

In addition, let if be a fixed bipartite graph on the same vertex set. In this paper we 
find precise formulae for the probabilities that H is edge-disjoint from G G B(s,t), that 
if is a subgraph of G, and that H is an induced subgraph of G. These probabilities are 
defined for the uniform distribution on B(s, t). In general, whenever we refer to a random 
element of a set, we always mean an element chosen uniformly at random. 

These formulae are obtained when the graphs in B(s, t) are sufficiently dense, the 
graph H is sufficiently sparse and the entries of s and t only vary within certain limits. 
The exact conditions are stated in Section [2j The starting point of the calculations is an 
enumeration of the set B(s, t, H) of graphs in B(s, t) which are edge-disjoint from H; see 
Theorem 12.11 

In the case m = n, the n x n binary matrix associated with the bipartite graph can 
also be interpreted as the adjacency matrix of a digraph which has no multiple edges but 
may have loops. By excluding the diagonal we obtain a parallel series of results for simple 
digraphs (digraphs without multiple edges or loops). These are presented in Section [3j 

These subgraph probabilities enable the development of a theory of random graphs 
and digraphs in these classes. As examples of computations made possible by this theory, 
we calculate the expected number of subgraphs isomorphic to a given regular subgraph. 
A particular case of interest is the permanent of a random 0-1 matrix with row sums s 
and column sums t. 

Now we briefly review the history of this problem. All previous precise asymp- 
totics were restricted to sparse graphs. Define g = max{s x , . . . , s m , ti, . . . , t n }, x = 
m£ix{x 1 , . . . , x m , jji, . . . , y n } and iV = ^\ Sj. Asymptotic estimates for bounded g were 
found by Bender [2] and Wormald [22]. This was extended by Bollobas and McKay [3] to 
the case g,x = 0(min{logm, logn} 1//3 ) and by McKay p2] to the case g 2 + xg = o(N 1 ^ 2 ). 
Estimates which are sometimes more widely applicable were given by McKay [llj . The 
best enumerative results for B(s,t) in the sparse domain appear in [8l [T7] . 

Although results about sparse digraphs with specified in-degree and out-degree se- 
quences can be deduced from the above, we are not aware of this having been done. Some 
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results using the pairings model have appeared [6]. For digraphs in the dense regime, some 
related work includes enumeration of tournaments by score sequence with possible forbid- 
den subgraph [131 EEl (HI E] , Eulerian digraphs [131 EE] > Eulerian oriented graphs [131 [21] , 
and digraphs with a given excess sequence [20] . 

For the case of dense bipartite graphs with specified degrees, an asymptotic formula for 
the case of empty H was given by Canfield and McKay [5] for semiregular graphs and by 
Canfield, Greenhill and McKay [3] for irregular graphs. The latter study is the inspiration 
for the present one. A similar study for graphs which are not necessarily bipartite is in 
preparation [13]. 

In related work using different methods, Barvinok [1] gives upper and lower bounds 
for \B(s,t, H)\ which hold very generally (from sparse to dense graphs) but which can 
differ by a factor of (mn) 0( * m+n \ Barvinok's results also give insight into the structure of a 
"typical" element of B(s, t, H), which he proves is close to a certain "maximum entropy" 
matrix. 

The paper is structured as follows. The results for bipartite graphs are presented in 
Section [2] and the corresponding results for digraphs can be found in Section [3j Then 
Section H] presents a proof of the fundamental enumeration result, Theorem 12.11 from 
which everything else follows. 

Throughout the paper, the asymptotic notation 0(f(m, n)) refers to the passage of m 
and n to oo. We also use a modified notation 0(f(m,n)), which is to be taken as a 
shorthand for 0(f(m,n)n°^ e ) , where the 0(1) factor is uniform over e provided e is 
small enough. 

2 Subgraphs of random bipartite graphs 

In this section we state our results for bipartite graphs. 

The starting point of the investigation is the enumeration formula given in the following 
theorem. Define m,n,s,t as in the Introduction and further define 



Note that s is the average degree on one side of the vertex bipartition, t is the average 
degree on the other side, and A is the edge density (the number of edges divided by mn) . 

Let H be a fixed bipartite graph on the same vertex set that defines B(s,t), namely 
{ui, . . . , u m } U {v±, . . . , v n }. For j = 1, . . . , m, k = 1, . . . , n, let Xj and y k be the degrees 



in 



n 
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of vertices Uj and v k of H, respectively, and further define 



5j = Sj - s + Xx j} r] k = t k -t + Xy k . 



Also define 



n 



j=l k=l jkeH 

m n 

j=l k=l 

In the case of Y and similar notation used in this section, the summation is over all 
j G {«!, . . . , u m } and k G {v 1: . . . , v n } such that UjV k is an edge of H. 

Theorem 2.1. For some e > 0, suppose that Sj — s, Xj, t k — t and y k are uniformly 
0{n l ' 2+e ) for 1 < j < m and 1 < k < n, and X = 0(n 1+2e ), for m, n — > oo. Let a, b > 
be constants such that a + b < |. Suppose that m, n —>■ oo with n = o{m 1+e ), m = o{n 1+e ) 
and 

1 - 2A) 2 / 5m 5n\ 

' H 1 < alogn. 



\ 6n 6m 
Then, provided e > is small enough, we have 

mn—X \ T-r / n—Xj \ -p-r / m—y k 



w-^i-CX^)"n(v)nC 

\ / j =1 \ i / k= i \ 

x expf-ifl - J?-) (i - C) _ -J— + (n 
\ 2\ 2AmnJ\ 2AmnJ 2Amn 



-b- 



The proof of Theorem 12.11 will be presented in Section HI As in the special case of 
empty H proved in p[], the formula for \B(s,t, H)\ has an intuitive interpretation. The 
first binomial and the two products of binomials are, respectively, the number of graphs 
with Xmn edges that avoid H, the number of such graphs with row sums s, and the number 
of such graphs with column sums t. Therefore, the exponential factor measures the non- 
independence of the events of having row sums s and having column sums t. Another 
expression for the product of binomials in the theorem is given below in equation (1TB]) . 



We can now employ Theorem EUJ to explore the uniform probability space over B(s,t) 
First we need a little more notation. For all nonnegative integers h, i define 

m n 

Rh,£ = 5^ S J X& V C h,£ = Y VkVk- 

j=l k=l 
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We will abbreviate R h0 = R h and C h0 = C h . Also note that R 1 = C 1 = XX and 
R Q1 = C 0A = X. Finally, let 

^1,1 — x jVk, = SjUki Yl,0 = x jVk- 

jk£H jkeH jk£H 

Theorem 2.2. Under the conditions of Theorem \2.1[ the following are true for a random 
graph G G B(s,t) provided e > is small enough: 

(i) the probability that G is edge-disjoint from H is (1 — X) x miss(m, to); 

(ii) the probability that G contains H as a subgraph is A x hit(rro, to), 

where 

XX /l In AX 2 1 C hl 



missira, to) = exp — — — | H — ; — — : — h 

1 2(l-A)Vn m/ 2(1 - A)mn 1-AVto m 

_ Y + A ^ , 2 + | A(1-2A) ^o, 3 + 

m 2 



X(l-X)mn 2(1 - A) V to m/ 6(1 - A) \ to 

1 — 2A ( R\,2 , ^1,2 \ 1 C*2,l 



2(1 — A)^ v to 2 m 2 / 2(1 — A) 2 \ n 2 m 2 



+ ^ -T— 7^ M+ + ^+ +0(TO 



and 



hit(m n) - exp A 1 ~ A ) X ("I + 1_\ + ^ ~ + I f ^ 1 . 1 + £i± 

\ 2A Vto m/ 2Xmn A V to rro 

_ -Lflhl + _ 1 + A , l + 2A /i?l, 2 C^2 

2A 2 U 2 m 2 i 2A Vto tot / 2A 2 V to 2 m 2 
(l + A)(l + 2A) ^o, 3 C , 3 \, ^-^0,1-^1,0 + ^1,1 , n , 

6A 2 U 2 + m 2 ) A(l — X)mn + U[U } 

Proof. The first probability in the statement of Theorem 12.21 is 

\B(s,t,H)\ 

\B(8,t)\ 

which can be expanded using Theorem 12.11 (One method is to apply ffl8l) below.) The 
second probability can be derived in similar fashion, or can be deduced from the first on 
noting that the probability that G includes H is the probability that the complement of 
G avoids H. □ 

In the standard model of random bipartite graphs on m + to vertices with expected 
edge density A, each of the mn possible edges is present independently with probability A. 
The probability that a random bipartite graph taken from the standard model is disjoint 



5 



from or contains a given set of X edges is (1 — X) x or X x , respectively. Therefore, the 
quantities miss(m, n) and hit(m, n) given in Theorem 12.21 can be interpreted as a measure 
of how far these probabilities differ in £>(s, t) compared to the standard model. Suppose 
that in addition to the conditions of Theorem 12. 2\ we also have 

X max \sj — s\ + XR 2 = o((l — A)n) , 

' ' ' (1) 

X max \t k —t\ + XC 2 = o((l — A)m) . 

k 

Then miss(m, n) = 1 + o(l). Similarly, if we have 

Imax |sj — s| + (1 — A)i2 2 = o(An), 

(2) 

Xmax |£ fc — £| + (1 — A)C 2 = o(Am). 

k 

then hit(m, n) = 1 + o(l). Requirements ([1]) and (jSJ) are both met, for example, if 
X = 0(n 1//2-2e ). Another interesting case is when s„- — s, Xj, t k — t and w fe are uniformly 
0(n e ) andX = 0(n 1 " 2£ ). 

To assist with the application of Theorem 12.21 we will give the simplifications that 
result when the graphs in B(s,t) are semiregular or when the graph H is semiregular. 

Corollary 2.1. In addition to the conditions of Theorem \2.Sl assume that s,- = s and 
t k = t for all j, k. Then 

XX ( 1 1 \ AX 2 XY 1 ! 

misslm, n) = exp ' 



1 1 N 
1 




n vtlj 


l + 2 


( R °> 2 + 




V n 


m 



X fR^ C^y A(2-A) ,Rw 0{n - b) \ 
-A)V n mJ 6 1-AfU 2 m 2 / J 



2(1- A) V n m / 6(1 - A) z 
and 

/ (1-A)X /1 lx (l-A)X 2 (1-A)y a)1 
mt(m, n) = exp 1 — I I + 



2 A \n m/ 2Xmn Xmn 



^& + ^)_l^& + ^) +0(n -»)). D 
2A V n m J 6A 2 V n 2 m 2 / / 



Corollary 2.2. In addition to the conditions of Theorem \2.Sl assume that Xj = x and 
Vk = V f or a H 3i k. (Note that Theorem \2.2\ requires x,y — 0(n 2e ) in that case.) Then 

/ X(xy-x-y) yR + xC Y _ b 

misslm, n) = exp ; — k -— -r h U(n ) 

1 2(1 - A) 2(1-A) 2 mn A(l - X)mn 

and 



\ f (1 - \){xy - x - y) yR + xC Y _ 6 , 

hitfm, n) = exp — — ■ — h U(n 

V 2A 2A 2 mn A(l - X)mn 

where Y = J2jkeii( s j ~ ~ *)■ D 



The next question we will address is the probability of H appearing as an induced 
subgraph. To be precise, suppose that H has no edges outside {u t , . . . , Uj} x {v x , . . . , v K } 
and let H JK denote the subgraph of H induced by those vertices. We will only consider 
the situation when the graphs in B(s,t) are semiregular. The corresponding result for 
irregular graphs can also be obtained using the same approach. 

The probability that H JK is an induced subgraph of G G B(s, t) is simpler to state in 
terms of some new variables. For I = 1, 2, 3, define 

J K 

j=l k=l 
Note that uii = u[ = X — \JK. 

Theorem 2.3. Adopt the assumptions of Theorem \2.1\ with Sj = s and t k = t for all j, k, 
and assume that J, K = 0{n l ^ 2+£ ). Then the probability that a random graph in B(s,t) 
has H JK as an induced subgraph is 

^jk-x ( (JK , (1-2A)wa/1 1\ u\ 



m n) AAmn 
{n + K)u 2 {m + J)u' 2 1 - 2A^ + ^3 \ + q^-* 



AAn 2 AAm 2 2AA 2 

Proof. Let H* be the complete bipartite graph on the parts . . . , Uj} and {v 1; . . . ,v K }. 
Then the probability that a random graph in B(s, t) has H JK as an induced subgraph is 

\B(s-x,t-y,H*)\ 
\B(s,t)\ 

This ratio can be estimated using Theorem 12.11 (or by combining Theorems 12. II and 12. 2\i . 

□ 

The argument of the exponential in Theorem 12.31 is o(l) if JK 2 = o(An) and J 2 K = 
o(Am). So, in those circumstances, the probabilities of induced subgraphs asymptotically 
match the standard bipartite random graph model for edge probability A. 

A related question asks for the distribution of the number of subgraphs of given type 
in a random graph in B(s,t). This deserves a serious study, which we will only just 
initiate here. A colour-preserving isomorphism of two bipartite graphs on . . . , u m } U 
{vi, . . . , v n } is an isomorphism that preserves the sets {ui, . . . , u m } and {v x , . . . , v n }. Let 
T(H) be the set of all graphs isomorphic to H by a colour-preserving isomorphism. We 
know that 



aut(if)' 
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where aut(H) is the number of colour-preserving automorphisms of H. 

When the graphs in B(s, t) are semiregular, the expected number of elements of T(H) 
that are contained in or edge-disjoint from a random graph in B(s, t) is clearly just \T(H) \ 
times the probability given by Theorem 12.21 and Corollary 12.11 

If this regularity condition does not hold, the calculation is more complex. Here we 
consider the case that the graph H is semiregular and leave the most general case for a 
future paper. 

We will need the following averaging lemma. 

Lemma 2.1. Let = {zf\ z£\ • • • , z$) be a vector in [—1, l] n such that Y^j=i z f^ = 0- 
Form z^\ z^ 2 \ ... as follows: for each r > 0, if z^ is the first of the smallest elements 
of and z[ is the first of the largest elements of z^\ then z^ +1 ^ is the same as z^ 
except that z^ r+1) and z[ +1) are both equal to (z^ + z { p)/2. Then z {n) E [-§, \} n . 

Proof. If z^ £ [— |, ^]™, then the fact that Y^j=i = implies that zf^ < and 
z[ > 0. Therefore z^ r+1 - ) has at least one fewer element outside [— |, |] than z^ does. 
The lemma follows. (In fact, ^LC 2 ™- 1 )/ 3 !) ^ [—h, |] n , but this improvement is not necessary 
for our application.) □ 

Theorem 2.4. Suppose that the conditions of Theorem \2.1\ apply with Xj = x and y k = y 
for all j, k. Then the following is true of a random graph G in B(s, t): 

(i) the expected number of graphs in 1(H) that are subgraphs of G is 

A \X(H)\exp[ — —2 + 0{n , 

V ^ 2A mn 

(ii) the expected number of graphs in T(H) that are edge-disjoint from G is 

fi \^iT/mi ( A(zy-x-y) yR + xC _ b 
(l-A) IK^Iexp^-^-^ — w - + (n ) 

Proof. Define z^°\ z^\ ... as in Lemma [2. 1\ with z^ = Sj — s for 1 < j < n. For r > 0, 
define 

m 

Y {r) (g, h) = Y^ z f T jth , where T j h = (t fch - t), 

and 

^ \ 2Amn 



j=l k: Uj v k £E(H) 
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For a permutation pair (g,h) G S m x S n , define H 9 ' h to be the isomorph of H with 
edge set {ufVjh \ UjV k G E(H)}. As (g, h) runs over S m x S n , each isomorph of H appears 
as H 9,h exactly ant(H) times. Therefore, by Corollary I2.2[ the expectation required in 
part (i) of the theorem is 

A X F (0) / (i-X)(xy-x-y) yR + xC _ b 
exp( — —2 h 0(n ) 



&ut(H) ^\ 2X 2Xmn 

For some r > 0, suppose that z^ r+1 ^ is formed from by averaging zf^ and zjf^ as 
in Lemma [2~T1 Then | g G S^} = if> n , so 



^-iE(-p(-^) + -p 



2Amn 



I Vex f ^m*> z / T J> \ f cx / 4 r)T i,fe + ^ r )T < 
2 ^ eXP V 2Amn M ex P^ 2Amn 



(g,h) 



eji 



+ exp (- 



z i 1 £,h + z l J-iM 



2 Aran 



V y^T (r+ljrp . {r+l)rp 

Y exp {- hm^l L 2± _ T «.* + *' T ^ + 0(n- 2 ^ 
^ V 2 Aran 2Aran 



(g,h) 



= F (r+1) exp(0(n- 2 )). 

By Lemma 12.11 there is some r = 0(n log n) such that z^ -' G [— n^ 1 ^ 2 , n~ l l 2 ] n . By the 
definition of F^ r °\ we have F 1 -'" ^ = ra\n\ exp(0(ra -1 )) , so F*- -* = ml nl exp(0(n -1 )) by 
induction. Part (i) of the theorem follows. Part (ii) is proved in identical fashion. □ 

A simple example of Theorem 12.41 at work is the enumeration of perfect matchings in 
the case m = n. Equivalently, this is the permanent of the corresponding n x n binary 
matrix. Most previous research has focussed on the case that the matrix has constant 
row and column sums. For s = t = o(ra 1 ^ 3 ), the asymptotic expectation and variance are 
known, while for s = t = n — 0(n 1_e ), the asymptotic expectation is known [3J. In the 
intermediate range of densities covered by the current paper, it appears that only bounds 
are known. The van der Waerden lower bound nl \ n (proved independently by Egorychev 
and Falikman) was improved by Gurvits [9] to s! f(s — l) s-1 /V~ 2 ) n s . The best upper 
bound is sl 1 ^ ~ nl A n+1 ^ 2A - ) (27rn) ( ' 1_A ^ ( ' 2A ' ) conjectured by Mine and proved by Bregman. 
See Timashev [IS] for references and discussion. 
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Applying Theorem I2.4( i) with x = y = 1 gives the following. 



Theorem 2.5. Suppose that m = n and s,t,X satisfy the requirements of Theorem \2.1[ 
Then the expected permanent of a random n x n matrix over {0, 1} with row sums s and 
column sums t is 

It is interesting to note that in the regular case R = C = 0, the average given in 
Theorem 12.51 is only higher than Gurvits' lower bound pj by a factor of A~ 1//2 (l + o(l)). 



3 Subdigraphs of random digraphs 

The adjacency matrix of a simple digraph is a square {0, l}-matrix with zero diagonal. 
Therefore, Theorem 12. II can be applied to enumerate digraphs with specified degrees, and 
the result can then be used to explore the corresponding uniform probability space. 

In this section, H denotes a fixed simple digraph on the vertices {w 1 , . . . ,w n }. Let 
V(s,t) be the set of all simple digraphs on vertices {w 1 , . . . , w n } with out-degrees s and 
in-degrees t, and let V(s, t, H) be the subset of T>(s, t) containing those digraphs which 
are arc-disjoint from H. 

For 1 < j < n let Xj, yj denote the out-degree and in-degree of vertex Wj in H, 
respectively. The quantities s — t, X, Sj, m, X, Y, and so forth are all defined by the 
same formulae as in Section [2] with m = n. In the definition of Y, the summation over 
jk G H should now be interpreted as summation over j,k G {1, . . . , n} such that WjW k is 
an arc of H. But note that A does not represent the arc-density of a digraph in T>(s,t). 
Instead the arc-density of a digraph in TJ(s, t) is given by 

p — s/(n — 1). 

We begin with the basic enumeration result for digraphs. 

Theorem 3.1. For some e > 0, suppose that s„- — s, Xj, tj — s and yj are uniformly 
O^n 1 /^ 6 ) for 1 < j, k < n, and X = 0{n 1+2e ), for n — > oo. Let a, b > be constants 
such that a + b < |. Suppose that n — > oo with 

(1-2A) 2 

< alogra. 



3A 
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Then, provided e > is small enough, we have 

\V(s,t,H)\ 

n 2 —X—n\ 1 j-T ( n—Xj—l\ /n—yj—1 
Xn * Jill s , [ t 



3=1 



3 / \ 3 



Proof. Let H be the bipartite graph obtained from H by replacing each vertex Wj by two 
vertices it,-, Vj, replacing each arc WjW k of H by the edge UjV k of H, and finally adding the 
perfect matching {ujVj \ j = 1, . . . , n} to the edge set of H. Then the degree sequences 
on the left and right of H and the total number of edges in H are given by 

Xj = Xj + 1, ifj = yj + 1, X = X + n, 

respectively (1 < j < n). The quantity Y for H satisfies 



Y = Y + ^2(s j -s)(t j -s) + 0{ 

3=1 



2-b\ 

n 



for any positive constant b < 1/2. Using this fact while applying Theorem 12.11 to H 
completes the proof. □ 

This formula for \T>(s, t, H) \ has an intuitive interpretation which is analogous to that 
given after Theorem 12.11 for the bipartite graph case. 

Using this enumeration theorem, we can explore the uniform probability space over 
T>(s,t). In each case, the proof is analogous to that of the corresponding theorem for 
bipartite graphs in Section [2j 

Theorem 3.2. Under the conditions of Theorem \3.1l the following are true for a random 
digraph G G T>(s, t) if e > is small enough: 

(i) the probability that G is arc- disjoint from H is (1 — p) x miss(n, n); 

(ii) the probability G contains H as a subdigraph is p X hit (n,n), 

where miss(m, n) and hit(m, n) are defined in Theorem \2.2[ □ 

The special cases of miss(m, n) and hit(m, n) provided by Corollaries 12.11 and 12.21 
apply here as well, as do the sufficient conditions ([I]) and ([2} for the probabilities in 
Theorem 13.21 to asymptotically match those in the standard random digraph model with 
arc probability p. 
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Next suppose that each arc of H has both ends in {w x , . . . , Wj}. Let Hj be the 
subdigraph of H induced by those vertices. For £=1,2, 3, define 



j=l k=l 

Note that Xi — Xi — X — pJ(J — 1). 

Theorem 3.3. Adopt the assumptions of Theorem \3.2\ with Sj = s and t k = t for all j, k, 
and assume that J = 0(n 1 ^ 2+£ ). The probability that a random digraph in T>(s, t) has Hj 
as an induced subdigraph is 

x/t \J(J-i)-x (J 2 , (!-2A)Xi Xi 
p (1 — p) v ; exp — 1 



2 



n 2An 4An 

(n + J)(x2 + X2) (l-2A)(x 3 + X3) , n( - b ,\ n 

o H Oln j . □ 

AAn 2 24 A 2 n 2 V ' J 

The argument of the exponential in Theorem l3.3l is o(l) if J 3 = o(An). So in that case, 
the probabilities of induced subdigraphs asymptotically matches the standard random 
digraph model for arc probability p. 

Let 1(H) be the isomorphism class of H and note that \I(H)\ = n!/aut(if), where 
aut(if ) is the number of automorphisms of H. By the same averaging technique as used 
to prove Theorem 12. 4[ we obtain the following. 

Theorem 3.4. Suppose that the conditions of Theorem \3.1\ apply with Xj = yj = x for 
all j . Then the following is true of a random digraph G in T>(s, t): 

(i) the expected number of digraphs in T(H) that are subgraphs of G is 

xwttv ( (1 - A) z(z - 2) (R + C)x _ 
p X \l(H)\ exp(-^ ^ >- - 2aV ; + 0(n 

(ii) the expected number of digraphs in T(H) that are arc- disjoint from G is 

xXi^/rrM ( \xix-1) (R + C)x _ b .\ _ 

(l- p )-|I W |exp(-^^__L_^ +0( „^. D 
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4 Proof of Theorem 



2.1 



In the remainder of the paper we give the proof of Theorem 12.11 The overall method 
and many of the calculations will parallel |4J, albeit with extra twists at each step, so we 
acknowledge our considerable debt to Rod Canfield. 

Outline of proof of Theorem \2.1\ The basic idea is to identify \B(s, t, H)\ as a coefficient 
in a multivariable generating function and to extract that coefficient using the saddle- 
point method. In Subsection 14. 1[ we write \B(s,t,H)\ = P(s,t, H)I(s,t, H), where 
P(s, t, H) is a rational expression and I(s, t, H) is an integral in m+n complex dimensions. 
Both depend on the location of the saddle point, which is the solution of some nonlinear 
equations. Those equations are solved in Subsection 14.21 and this leads to the value of 
P(s,t,H) in ( fl9l) . In Subsections I4.31 - f4.6l the integral I(s,t,H) is estimated in a small 
region 7Z' defined in ([30]) . The result is given by Lemma [4.31 together with ([22]) . Finally, 
in Subsection 14.71 it is shown that the integral I(s,t,H) restricted to the exterior of 7Z' 
is negligible. Theorem [O then follows from (jl), ([19]), Lemmas OHO and ([22]). □ 

We will use a shorthand notation for summation over doubly subscripted variables. 
If Zj k is a variable for 1 < j < m and 1 < k < n, then 



rn 



• ^ ^ Zjki ^»k ^ ^ Z jk i ^ ^ ^ ^ ^jki 

k=l j=l j=l k=l 

n— 1 m— 1 m— 1 n— 1 

^j* ^ , Zjki ^*k ^ > ^jki ^ , ^ , ^jki 

k=l j=l j=l k=l 

for 1 < j < m and 1 < k < n. 

For 1 < j < m and 1 < k < n define hj k = 1 if UjV k is an edge of H and hj k = 
otherwise. Then define the sets 

Xj = { k | 1 < k < n, h jk = 1 }, Xj , = { k \ 1 < k < n, h jk = }, 
Y k = { j | 1 < j < n, h jk = 1 }, Y k = { j | 1 < j < n, h jk = }, 

for 1 < j < m and 1 < k < n. 

The notations J2jkeH anc ^ J2jkeH indicate sums over the sets {(j,k) \ 1 < j < m, 
1 < k < n,hj k = 1} and {(j, k) \ 1 < j < m,l < k < n,hj k = 0}, respectively, and 



13 



similarly for products. We also define summations whose domain is limited by H. 

Z 3»\H = Z 3 k '> Z »k\H — z jki z »m\H = z jki 

kex 3 jeY k // 

Z j*\B = Z i kl Z »k\H = Z i k ' Z »»\H = Z 3 k ' 

hex, jer k jken 



Under the assumptions of Theorem 12.11 we have m = 0(n) and n = 0{m). We also 



Cl^C 2 +C 3 £ C 4 +C 5 £ 



n ' 



have that 8 < A 1 < 0(\ogn), so A 1 = 0(1). More generally, A Cl m 
0(n C2+Ci ) if c^c^c^c^c^ are constants. 

We now show that the assumptions of Theorem 12.11 imply that 

m = o(A 2 n 1+£ ), n = o(A 2 m 1+£ ). (3) 

If A > ^ then ([3D follows immediately. If A < ^ then (1 - 2A) 2 > \ and so the 
assumptions of Theorem 12.11 imply that 1/A = 0(\ogn/{m/n + n/m)) . This implies (EJ). 

4.1 Expressing the desired quantity as an integral 

In this section we express \B(s, t, H) \ as a contour integral in (m + n)-dimensional complex 
space, then begin to estimate its value using the saddle-point method. 

Firstly, notice that \B(s, t, H) | is the coefficient of ul 1 ■ ■ ■ w^ 1 ■ ■ ■ in the function 

Y[ (1 + UjWk). 
jkGH 

By Cauchy's coefficient theorem this equals 



\B(s,t,H)\ = - (p ■■■ & — — - — 3 e J — T du 1 ---du m dw 1 ---dw n , 



(2m) m+n u^ +1 ■ ■ ■ u s - +1 w^ +1 ■ ■ ■ wy 1 

where each integral is along a simple closed contour enclosing the origin anticlockwise. 
It will suffice to take each contour to be a circle; specifically, we will write 

Uj = qje 3 and w k = r k e Vk 

for 1 < j < m and 1 < k < n. Also define 

A jk 



1 + qf k 
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for 1 < j < m and 1 < k < n. Then \B(s, t,H)\ = P(s, t, H)I(s, t, H) where 



P(s,t,H) = - , , , . 

(4) 



r r n,^ff(i+A 1fe (e i( ^ fc) -i)) 

J — 7t ./-* exp(z ^ i=1 s^- + z )^ fc=1 t fc fc ) 

= (6*!, . . . , 9 m ) and <j> = (<p 1} ...,<p n ). 

We will choose the radii qj, r k so that there is no linear term in the logarithm of the 
integrand of I(s, t, H) when expanded for small 6, <f>. This gives the equation 

m n 
jkeH J=l k=l 

For this to hold for all 6, 0, we require 

X j.\R = s j (1 < J < m), 



Kk\R = h {l<k< n). 



(5) 



The quantities Xj k have an interesting interpretation. If edge UjV k is chosen with 
probability A jfc independently for all j, k E H, then the expected degrees are s, t. 

In addition to the quantities defined before the statement of Theorem 12.21 we define 
for j = 1, . . . , m, k = 1, ... ,n, 



4.2 Locating the saddle-point 

In this subsection we solve (jSD and derive some of the consequences of the solution. As 
with the whole paper, we work under the assumptions of Theorem 12.11 

Change variables to {a,j}™ = i, {b k } k=1 as follows: 

1 + dj l + b k 

Qj = r ~, — > r k = r ~ 2T> (6) 

1 — r a j 1 — r b k 

where 

A 



1 - A 

Equation (JSJ) is slightly underdetermined, which we will exploit to impose an addi- 
tional condition. If {qj}, {r k } satisfy (jSJ) and c > is a constant, then {cqj}, {r k /c} also 
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satisfy (JSJ). From this we can see that, if there is a solution to (jSJ) at all, there is one 
for which X^=i( n — x j) a j < an d X]fc=i( m ~~ Uk)^k > 0> an d a l so a solution for which 
Sj=i( n — x i) a i > an d Sfc=i( m ~ Z/fc)^ < 0. It follows from the Intermediate Value 
Theorem that there is a solution for which 

m n 

"^(n - x 3 )dj = J^(m - y k )b k , (7) 

3=1 k=l 

so we will seek a common solution to (jSJ) and (171) . 
From (jSJ) we find that 

Ajjfe/ A = 1 + dj + b k + Z jk , (8) 

where 

a,-fei.(l — r — r a~ — r b k ) 

Zjk = Jk{ - ~ 2 - J (9) 

1 + r ajOfc 

and that equations flS} can be rewritten as 



-l = (n- Xj )aj + h + 



Am 2mn m 



mn m m 2mn 
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(10) 



Summing (ITU1) over all j, k, respectively, we find in both cases that that 

m n 

X = z^2( n ~ x i) a i + 5Z( m ~~ yk "> bk + Z -\B ■ ( n ) 

3=1 k=l 

Equations <£7§ and ffTTl) together imply that 

m n 

^2(n - Xj.)^. = J^(m - y k )b k = \{X - Z,^ R ) . 

3=1 k=l 

Substituting back into (flUl) . we obtain 

a j = &j\ a lj • • • j a m; &1> ■ ■ ■ j b n ), 

b k = B fc (a 1 , a m , 6„), 

for 1 < j < m, 1 < k < n, where 

. . dj X &jXj 

Aj a 1; . . . , a m , 6 n = - + 

An 2mn n 

_ ELi Vkh ^keXj b k _ Zj, lR Z„ lB 

mn n n 2mn 



(12) 



By the same argument as in [3], equation (T12I) defines a convergent iteration starting 
with a,j = b k = for all j, k. Four iterations give the following estimate of aj. The value 
of b k follows by symmetry, while Zj k follows from (Q. 

dj dj.i-j dj.rj dj.r'j X (1 - 2X)5 j X d 2 X 



An An 2 An 3 An 4 2mn AAmn AAmn 3 
XjX x)X A(7-10A)X 2 3(1 - 2X)5 jXj X (1 - 2X)Y 

mn mn 16Am n AAmn 4Av4m n 

<LC 2 (1-2A)5 2 C 2 (Lx-Ca 3XC 2 X# 2 



_i_ _; ^ — i _i_ 

. . . 



2XAm n AXA' ''m 2 n 3 XAm 2 n 3 8Am 3 n 2 8Am 2 n 3 
(1-2X)R 2 C 2 XjR hl C 1A (l-2X)5jC hl Xj C hl 



8AA 2 m 3 n 3 Amn 3 Xm 2 n 2XAm 2 n 2 Xm 2 n 2 
Y 0tl C lj2 _Jj_ (1 - 2X)6 j J j _ ^J,- ./>/,• 



Xm 2 n 2 Am 3 n Amn 2AAmn 2 2AAmn 3 Amn 2 
XjJ 3 - (1 - 2X)5 j x j J j 3(1 - 2A)XJ i (1 - 2A) J 2 



j ./ J j _^ ' J J J v ^ J_ j 

Amn 3 AAmn 3 4Am 2 n 2 2XAm 2 n 2 



( -R 2 C 2 \ Jj 1 r- 1 V ^ 2 

( — + — ) w 2 2 + E <V%' + E ™ 

v n m / 2XAm n Am n Am n ~~ 



X ^ J j 2 (1-2A) 



Ef*" oulu E ^ + » E <V^' 



m 2 n 2 "™ 2AAm z n z "" 2AAm z n , , 

keXi kex, (/ ifc ') 



1 — 2A)<5j Xj 



+ ( v^rr + , 2 3 + j ( n 2^ ^ + m E <V ) 

Mm n 2XAm n Am n / \ t - ; j / 

+ iE <W + t-t-2 E E v + 5 (^" 5/2 )' 



Amn "7^ Am n , , „ 

(j',k') {j',k')k"ex., 



where the notation E(/,fc') means Efc'ex, E/eK/- 

A sufficient approximation of Xj k is given by substituting this estimate into (jHJ). In eval- 
uating the integral I(s,t,H), the following approximations will be required: 



A , (1 _ A , ) . A(1 _ A)+ (iz^ + (iz|^_4_4 

n m n m 

(l-12A)5 jVk (l-2X)S j x j (l-2X) Vk y k 



2Amn n 2 m 2 



+ (l-2X)(J j + K k -XX) + ~ 3/2) 
mn 
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A ifc (l - X jk )(l - 2X jk ) = A(l - A)(l - 2A) + 



(1 - 12^)^- 



n 



(14) 



+ v L!ll + o(n 1 , 

m 

A ifc (l - X jk ){l - 6X jk + 6X%) = A(l - A)(l - 12A) + d( n - 1/2 ). (15) 
We now estimate the factor P(s,t, H). If 



a = n ^(i-^) 1 - 

jk&H 



then 



a- = n P^TW^' 

n nr w) 

;;■-// .7 = 1 fc=l 



jkeH 

m n 

= n (i+^n^n^ 

using (jSJ). Therefore, the factor P(s,t, H) in front of the integral in (j3J) is given by 

P(s,t,H) = (2n)- {rn+n) A' 1 . 
We proceed to estimate A. Writing Xj k = A(l + we have 



x\i-xy- x j J ^ °vi-a 

A ^ 2 A(l - 2A) _ 3 , A(l - 3A + 3A 2 ) _ 4 , ^ f z% 



^~ " ' * ^ \ Z jk , \\2 Z jk~^~ 10/1 \ N 3 ^ 



2(1 -A) ^ 6(1 -A) 2 3 12(1 - A) J V(l-A) 



(16) 



We know from (jSJ) that A..|# = Amn, which implies that = X, hence the first term 
on the right side of fflBI) contributes A AX (1 — A)~ AX to A. Now using (jSJ) we can write 
ZjT- = cij + b k + Zj k and apply the above estimates to obtain 



A= (A A (1-A) 1 - A ) mn (l-A)- A 



, R 2 C 2 R 2 C 2 X l X l (1-2X)/R 3 C 3 \ 

xexp — - + —. — + — 2 1 \^ + —^) 17 

• \AAn AAm 8A 2 m 2 n 2 4 Amn 24A 2 V n 2 m 2 J 

(1 - 6A) (R 1+ C ± \ + _J_ + J^L + + Ofn- 1 / 2 )^) 

96 A 3 W 3 m 3 / 2Amn 4An 2 4Am 2 / 
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As in [4 J , our answer will be simpler when written in terms of binomial coefficients. 
Using an accurate approximation of the binomial coefficients (such as [U Equation 18]), 
we obtain that 

mn-XY 1 yr ( n-xA A / m-y k \ (A A (1 - A) 1 ~ A )~ mn (l - X) x 

Xmn J 11 ^ Sj J 11 V h J (47rA) (m + n-l)/2 m (n-l)/2 n (m-l)/2 

f R 2 C 2 \-2Arm n\ 1-4A/R 2 C 2 
6XP \ 4:An AAm 24 A \n m) 16 A 2 \n 2 m 

2AA 2 \n 2 m 2 ) 96A 3 V ra 3 m 3 
A 2 X(m + n + X) _ -R 2 ,i _ C 2A Q( n ~y* 
iAmn AAn AAm 2 

Putting (TTTI) and (ITS]) together, we find that 
P(s,t, H) = k-\2it)- {m+n) 

A (m+n-i)/2 (n-l)/2 (m-l)/2 / vr\-l m / \ n / 

A y " nv " rv " (mn—X\ -r-r (n—Xj \ -i-r lm—y k 



27r"" + " +1,/2 V h»» J HI s* J n I /, 



j=l v "J 7 fe=l V ' fe 

, 1 - -M //,/ | n\ # 2 C 2 l-4A /R 2 | C 2 
X<X1> 1 24A \n m) 8A 2 m 2 n 2 16 A 2 V n 2 m 2 



A A \m n) 2 Aran 



(19) 



4.3 Evaluating the integral 

Our next task is to evaluate the integral I(s, t, H) given (J4]). 

Let C be the ring of real numbers modulo 2tc, which we can interpret as points on a 
circle, and let z be the canonical mapping from C to the real interval (— n,7r]. An open 
half-circle is C t = (t — tt/2, t + vr/2) C C for some £. Now define 

C N = { v = (v 1 , . . . , v N ) G \v 1 ,...,v N E C t for some teR}. 

If v = (vi, . . . , v N ) G C(f then define 

More generally, if v G C t then define v = t + (t> 1 — t, . . . , v N — t). The function v — > v 
is well-defined and continuous for v G C^. 
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Let TZ denote the set of vector pairs (0, <fi) £ C m x C n such that 

\e + 0| < (ffin)- 1/2+2£ , 

|^-|<n~ 1/2+e (l<j<m), (20) 
|4l < m~ 1/2+£ (1 < Jfe < n), 



where 9a = 9 a — and <\> k = (f) k — (f>. In this definition, values are considered in C. The 
constant e is the sufficiently-small value required by Theorem 12.11 

Let I n "(s,t, H) denote the integral I(s,t,H) restricted to any region TZ . In this 
subsection, we estimate I n '(s,t, H) in a certain region 71' D 71. In Subsection O wc 
will show that the remaining parts of I(s,t,H) are negligible. We begin by analysing 
the integrand in TZ, but for future use when we expand the region to 7Z (to be defined 
in ( |30j) ). note that all the approximations we establish for the integrand in TZ also hold in 
the superset of TZ' defined by 

\0 + 4>\ < 3(mn)- 1/2+2£ , 

\0j\ < 3n~ 1/2+£ (l<j<m-l), 

\L\ < 2n- 1 / 2+3£ , (21) 
\4> k \ < 3m~ 1/2+£ (l<k<n-l), 
|AJ < 2m~^ £ . 

Define = (9 1 , . . . ,9 m _i) and (f) = (0 1; . . . , n _i). Let T x be the transformation 
Ti(0, <j>, v, ■0) = (6, cf>) defined by 

v = + (j), i/j = 0-<j>, 

together with 0j = 9j — (1 < j < m — 1) and (f) k = (j) k — <f> (1 < k < n — 1). We also 
define the 1-many transformation T* by 

Tl{0,&v) = {jT x {b,&v^). 

After applying the transformation T x to In(s, t, H), the new integrand is easily seen to 
be independent of i/j, so we can multiply by the range of i/j and remove it as an independent 
variable. Therefore, we can continue with an (m + n — l)-dimensional integral over S such 
that TZ = T?(S). More generally, if S" C (-|vr, |vr) m+n " 2 x (-27r,2vr] and TZ" = T?(S"), 
we have 

I n „ (s, t, H) = 2vrmn / G{0, <j>, v) d0d(j>du, (22) 
Js" 
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where G(0, 0, v) = F(7i(0, 0, »/, 0)) with F(0, 0) being the integrand of I(s, t, H). The 
factor 27imn combines the range of ip, which is An, and the Jacobian of T 1; which is mn/2. 

Note that 5 is defined by the same inequalities ( l20i) as define 7?-. The first inequality 
is now \u\ < (mn) _1 ' 2+2E and the bounds on 



m— 1 n— 1 

0fc 



#m = - X] ^' alld ^ m = ~ E 

j=l k=l 



still apply even though these are no longer variables of integration. 
In the region S, the integrand of ([221) can be expanded as 



y 3 



G(0, fa v) = exp ( - ^ (A + a jk ){v + 9^ + 4> k f - i ( A 3 + /%fc)0 + h 

jk<=H jkeH 

+ ^(A 4 + 7 , fe )(^ + ^ + 4) 4 + o(A W + o 3 + k\ 5 )) 

jkeR jkeH 

(m n m n 
- Y Y( A + a 3^ v + h + - * Y E^ + ^-fc)^ + ^ + 

3=1 fc=l jr = l fc=l 

m n \ 

+ E Y^ + + ^ + + E A ^ + § i + + o(^^ 1/2 ) ) . 

j=i fc=i jkeH ' 

Here o; jfc , and jj k are defined by 

|^ife(l ~~ ^ife) = ^4 + ajfc; 
|A jfc (l - A ifc )(l - 2A ife ) =A 3 + /3 jk , (23) 
^A jfc (l — A jfc )(l — 6A jfc + 6A^ fc ) = A 4 + 7 jfc , 

where 

A = lA(l - A), A 3 = iA(l - A)(l - 2A), and A 4 = ^A(l - A)(l - 6A + 6A 2 ). 

Approximations for /J^, 7^ were given in (|T51) - (TT51) . Note that aj k in this paper is 
slightly different from in [I], but it is still true that atj k , Pjkiljk = 0{n~ 1 ^ 2 ) uniformly 
over j, k. 
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4.4 Another change of variables 

We now make a second change of variables (0, 0, z/) — ^(C? where £ — (d, . . . , Cm-i) 
and £ = (£i, . . . , £ n _i), whose purpose is to almost diagonalize the quadratic part of G. 
The diagonalization will be completed in the next subsection. The transformation T 2 is 
defined as follows. For 1 < j < m — 1 and 1 < k < n — 1 let 

h = Cj + ctti, <Pk = £k + dpx, 

where 

c = ^ and d = m 

m + m ' n + n 

and, for 1 < h < 4, 

m— 1 n— 1 

j=l k=l 

The Jacobian of the transformation is (mn) -1 ' 2 . In [5J, this transformation was seen 
to exactly diagonalize the quadratic part of the integrand in the semiregular case. In 
the present irregular case, the diagonalization is no longer exact but still provides useful 
progress. 

By summing the equations 9j = (j + C7r i an d 4>k — £k + dpi, we find that 

m— 1 

3=1 
n-1 

Pi = n 1/2 ^2<f) k , 
k=l 

where the inequalities come from the bounds on 8 m and n . This implies that 

^■ = ^■ + 0(0 (l<j<m-l), 
^ k = 4> k + d(n- 1 ) (l<k<n-l). 

The transformed region of integration is T 2 ~ 1 (iS), but for convenience we will expand it a 
little to be the region defined by the inequalities 

| C .| < | n -V2+e (l<j< m -l), 

\ik\<\m l/2+£ (l<*<n-l), 

In^Km^n- 1 ^, (25) 

I I - 1/2 -l/2+£ 
11^/ \-1/2+2e 

|z/| < (mn) ' 



I I - 1/2 -l/2+e 



|pi| < n 1 ' 2 m 1 / 2+£ ; 



(24) 
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We now consider the new integrand E 1 = exp(L 1 ) = G o T 2 . As in [5], the semiregular 
parts of L l (those not involving <x,- fc , (3j k , 7,- fe or H) transform to 

— Amnv 2 — Ann 2 — Amp 2 — 3%A 3 nvn 2 — 3iA 3 mvp 2 + GA A ir 2 p 2 

— iA 3 nir 3 — iA 3 np 3 — 3iA 3 cn7i 1 7i 2 — 3iA 3 dmp l p 2 + A A nn A + A A mp A + 0{rT 1 ^ 2 ). 

To see the effect of the transformation on the irregular parts of the integrand, write 
Cm = 6 m — C7r i an d in = @n — dpi- From (12^1) we can see that C, m = 0{rT 1 ^ 2 ) and 
£ n = 0{rT l l 2 ). Thus we have, for all 1 < j < m and 1 < k < n, (j + £ fc = 0{rT 1 ^ 2 ) and 
cn 1 + dp 1 + v = 0(rT l ). Recalling also that otj k , /5L- fc , jj k = 0(n -1 ^ 2 ), we have 



EEm^^ 

j=l k=l m n 

= E E a A(Cj + + 2 (Ci + Zk)(v + cvri + dp x )) + 0(n- 1/2 ), 

3=1 k=l 



EE^ + ^ + = EEw + + o( 

j=l k=l j=l k=l 



n 1/2 ) 



m n 



j=i k=i 



n 1/2 ) 



Moreover, the terms on the right sides of the above that involve ( m or £ n contribute only 
0{n~ 1 ^ 2 ) in total, so we can drop them. Combining this with (I26I) . we have 

L± = —Amnv — Anix 2 — Amp 2 — 3iA 3 nvir 2 — 3iA 3 mvp 2 + 6A A ir 2 p 2 

— iA 3 riTr 3 — iA 3 np 3 — 3iA 3 cmi 1 7i 2 — 3iA 3 dmp 1 p 2 + A A nir A + A A mp A 

m— 1 n—1 

-EE a i*(«i + + + + c7r i + ^i)) ( 27 ) 
i=i k=i 

m—l n—1 

- » E E + + A E «j + + o(^ 1/2 ). 

j=i fc=i jfeeif 
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4.5 Completing the diagonalization 

The quadratic form in L 1 is the following function of the m + n — 1 variables £, v: 
Q = —Aranv 2 — An7r 2 — Amp 2 + A ((^ + £ fc ) 2 

m-1 n-1 (28) 

M(C; + + 2(0 + Cfc)(^ + cvtx + 

i=i fc=i 

We will make a third change of variables, (£, v) = T 3 (a,r,fi), that diagonalizes this 
quadratic form, where er = (a 1 , . . . , er m _i) and r = (r 1; . . . , r n _ 1 ). This is achieved using 
a slight extension of [T6l Lemma 3.2]. 

Lemma 4.1. Let U and Y be square matrices of the same order, such that U~ l exists 
and all the eigenvalues of U~ X Y are less than 1 in absolute value. Then 

(i + Yu~ l y 112 (u + y) (i + u- x y)- 1/2 = U, 

where the fractional powers are defined by the binomial expansion. □ 

Note that U~ 1 Y and YU^ 1 have the same eigenvalues, so the eigenvalue condition 
on U~ 1 Y applies equally to YU^ 1 . If we also have that both U and Y are symmetric, 
then (I + YU^ 1 )^ 1 ^ 2 is the transpose of (J + U^Yy 1 ^ 2 , as proved in [1]. Let V be the 
symmetric matrix associated with the quadratic form Q. Write V = V d + V nd where 
has all off-diagonal entries equal to zero and matches V on the diagonal entries, and V nd 
has all diagonal entries zero and matches V on the off-diagonal entries. We will apply 
Lemma [4. II with U = V d and Y = V nd . Note that V d is invertible and that both V d and 
Vn d are symmetric. Let T 3 be the transformation given by T 3 (cr,T,fi) T = {C,£,v) T = 
(/ + V d ~ 1 V nd )~ 1//2 (cr, t, /x) T . If the eigenvalue condition of Lemma H~T1 is satisfied then this 
transformation diagonalizes the quadratic form Q, keeping the diagonal entries unchanged. 

From the formula for Q we extract the following coefficients, which tell us the diagonal 
and off-diagonal entries of V. Define x'j = Xj — hj n for 1 < j < m — 1, and y'k = Vk~ h-mk 
for 1 < k < n — 1. Then: 

[C|] Q = -An - (1 + 2c)a i# + Ax'j, 
[£} Q = -Am - (1 + 2d)a, k + Ay' k , 
[u 2 ] Q = —Amn, 
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[C'iCjJ Q = -2c( a ii* + «; 2 *) til ^ J2), 

[Cj£k] Q = ~ 2a jk - 2da^ - 2ca* k + 2Ah jk , 
fe^fcj Q = -2rf(a* fcl + a* k2 ) (fci ^ fc 2 ), 

Using these equations we find that all off-diagonal entries of V^f 1 l^ ld are 0(n~ 3//2 ), except 
for the column corresponding to u, which has off-diagonal entries of size 0{n ), and 
the entries corresponding to for /i^ = 1, which have size 0{n~ l ). Similarly, the 
off-diagonal entries of V^ d V^ _1 are all 0(n~ 3 ^ 2 ), except for the row corresponding to v, 
which has off-diagonal entries of size 0{n~ 1 ^ 2 ), and the entries corresponding to QCk f° r 
hj k = 1, which have size 0{pT l ). To see that these conditions imply that the eigenvalues 
of V d _1 V nd are less than one, recall that the value of any matrix norm is greater than or 
equal to the greatest absolute value of an eigenvalue. The oo-norm (maximum row sum 
of absolute values) of V d _1 V nd is 0{n~ 1 ^ 2 ) 1 so the eigenvalues are all 0(n -1 ^ 2 ). 
We also need to know the Jacobian of the transformation T 3 . 

Lemma 4.2 (@j). Let M be a matrix of order 0(m + n) with all eigenvalues uniformly 
0(n -1 ^ 2 ). Then 

det(J + M) = exp(trM-±trM 2 + 0(n~ 1/2 )). □ 

Let M = Kf 1 ^ / Qd- As noted before, the eigenvalues of M are all Ofn -1 / 2 ) so 
Lemma 14.21 applies. Noting that tr(JVf) = and calculating that tr(iW 2 ) = 0(n~ l ), 
we conclude that the Jacobian of T 3 is 

det ((/ + MY 1 ' 2 ) = (det(J + M)) ~ 1/2 = 1 + 0(r^ 1/2 ). 

To derive T 3 explicitly, we can expand (I + V d ~ 1 V^ ld ) -1 ^ 2 while noting that ctj* = 
0(n 1 ' 2+e ) for all j, = 0{m l l 2+e ) for all k, = 0{mn 2e + nm 2e ), R < mn 1+2e and 

n ^ l+2e 

This gives 

^ = ^ + 22{ 2An + ° {n ] h' 

^ /a jk + daj* + ca* k ~ _ 2 n\ c . f a j* . PC/ . PC/ -3/2\ 

fc=i 
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m—1 



= & + ^ +0(n ))Q 

+ E + 5 («- 2 ))4' + + 8(„-'))„ + 6 (B -*»), 

fc'=l 

m—1 Ti—l 
3=1 k=l 

for l<j<m — 1, 1 < /c < n — 1. 

The transformation T 3 _1 perturbs the region of integration in an irregular fashion that 
we must bound. From the explicit form of T 3 above, we have 



m—1 n— 1 



°i = Cj + E °( n ~ V %' + E °^ 12 )^ + 0(n^/2 )zy + 0(n -s/ a) = g + 0(n -i ); 

/=! fc=l 
m—1 n— 1 

r fe = £ fe + E 5 (^~ 3/2 )C, + E 0^' 2 )Z k > + 0(n- 1/2 )z/ + 0(n- 3/2 ) = 4 + O^ 1 ) 
i=i fc'=i 

for 1 < j < m — 1, 1 < k < n — 1, so <r, r are only slightly different from £, 
For // versus z/ we have 

fi = v + 0(7T 1+2£ A4) + 0(m^ 1+2 7A) 
= u + o((mra) _1 ^ 2+2e ) , 

where the second step requires ([3]). This shows that the bound \v\ < (mn) _1 ' 2+2£ is 
adequately covered by |//| < 2(mn)~ 1//2+2e . 
For 1 < /i < 4, define 

m—1 n— 1 

/^/ l = E cr ^' ^ = E Tfe/l - 

j=i *;=i 

From (I25j) . we see that [vri | < m 1/,2 n _1 ^ 2+£ and |px| < m _1 ^ 2+e n 1//2 are the remaining 
constraints that define the region of integration. We next apply these constraints to 
bound and v x . From the explicit form of T 3 , we have 

m—1 m—1 / i 

'C a,-* + a- 



^ = ^ + EE( 2Ara )JC/ 

m—1 n— 1 j m—1 

E E C Jt 2 1 — + 5 <' r2 >)«* + E + si"- 1 ))- + d(n-^ 



3=1 k=l j=l 
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COI** 1/2 -l/2+e . -1/2+e 1/2 . a ** 

= 7Ti H — m n H ; — m ' n ' H — v 

1 2An 2An 2An 

n—l , \ m— 1 

+ (' + *» " D) E ^ + ^ E <V.C, + 

fc=l /=! 

= 'i + £ W i E«/.C/ + 0(n-^) (29) 
/=i 

= ttj + 0(A _1 mn -1+2e ) 

/ 1/2 -l/2+5e/2\ 

= 7Ti + o(m ; n ' ' J. 

To derive the above we have used 1 + c(m — 1) = m 1//2 and the bounds we have established 
on the various variables. For the last step, we need ([3]), which implies that A" 1 mn~~ 1+2e = 
oim^n- 1 ' 2 ^ 2 ). 

Since our region of integration has < m 1 ^ 2 n~ 1 ^ 2+£ , we see that this implies the 
bound < m 1//2 n~ 1//2+3e . By a parallel argument, we have 

, / -l/2+5e/2 1/2 \ 
V X — Pi + °( m n J) 

which implies |^i| < n 1//2 m~ 1 ^ 2+3£ . Putting together all the bounds we have derived, we 
see that 

T.i\T 2 \S))CQnM, 

where 

Q = { | Oil < 2n~ 1/2+£ , j = 1, . . . , m - 1 } n { \r k \ < 2m~ 1/2+£ , k = 1, . . . , n - 1 } 
n{M < 2(mn)- 1/2+2e }, 
M = { |^| < m 1/2 n- 1/2+3£ } D { hi < n 1/2 m- 1/2+3£ }. 

Now define 

s' = T 2 (T 3 (gn.M)), 

K' = IKS'). 

We have proved that 5 3 5. Also notice that 1Z is contained in the region defined by 
the inequalities (12T|) . As we forecast at that time, our estimates of the integrand have 
been valid inside this expanded region. It remains to apply the transformation T 3 _1 to the 
integrand (1271) so that we have it in terms of (cr, r, /i). The explicit form of T 3 _1 is similar 
to the explicit form for T 3 , namely: 

m— 1 / i \ n—l , 

'C(a,-„ + a,' ) ~ _ 9 \ ^ /Qtwi. + dot,-* + ca» 



-(l£ +6 <"~ , >> +6 <"- 3/2 >- 
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& = * - e ( a,t "S; +ca,t + 6 ^ 2 >H- - e + 5 <"- 2 )) v 
' -(^+ s <»>+ 5 <»-V J 

m— 1 n—1 
3=1 k=l 

for 1 < j < in— I, 1 < k < n — 1. In addition to the relationships between the old and new 
variables that we proved before, we can note that n 2 = fi 2 + 0(n~ 1 ^ 2 ), p 2 = v 2 + 0(n -1 ^ 2 ), 
7r 3 = /i 3 + 0(n -1 ), p 3 = z/ 3 + 0(n~ l ), 7r 4 = /i 4 + 0(n~ 3//2 ), and p 4 = z/ 4 + 0(n~ 3//2 ). 

The quadratic part of L 1; which we called Q in ( ]28l . loses its off-diagonal parts ac- 
cording to our design of T 3 . Thus, what remains is 

to— 1 n—1 

—Amnp 2 — (An + (1 + 2c) Oj* — Ax'j)a 2 — ^(Am + (1 + 2d)a* k — Ay'^rl 

3=1 k=l 

= —Amnp 2 — Anp 2 — Amv 2 

to— 1 n—1 

- - Ar^o-J - ^(a* fc - Ay' k )rk + 0{n~ 1/2 ). 

3=1 k=l 

Next consider the cubic terms of L\. These are 

— 3iA 3 nvi\ 2 — 3%A 3 mvp 2 — iA 3 mc 3 — iA 3 np 3 

to— 1 n—1 

- 3iA 3 cnn 1 'K 2 - SiA^np^ - * E X] + & 

3=1 k=l 

We calculate the following in Q(l A4: 

• to— 1 n—1 

-1/2^ 



,fc) 3 - 



3 '^4. / TO— 1 7t— -L \ 

-3zA 3 raz/7r 2 = -3iA 3 n/iu 2 + 3 2 f o^cx,- + a^r*. j + 0( 

^ i=i fc=i ' 

. . 3zA 3 / k 2 

-2A 3 n7r 3 = -«A 3 n/i 3 + ^ 2^ c(c-* + ayj^ay, 



m— 1 n—1 



//* — ± n — _l \ 

+ E E( a J fc + + ca *k) a ] T k J + <5(rT 1/2 ), 

7=1 k=l ' 



O A 2 m-1 

6iA 3 c mp 2 \ ~ _i/2 N 

-ZiA z cmv 1 'K 2 = -3%A 3 cnp x p 2 H — + °( n )> ( 31 ) 



TO— 1 n—1 to— 1 n—1 

E E + ^) 3 = -» E E M°i + r *) 3 + o{n- l/2 ), 

j=l k=l 3=1 k=l 
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and the remaining cubic terms are each parallel to one of those. The proof of (I3T]) is 
similar to the proof of ( |29i) . 

Finally we come to the quartic part of L 1 , which is 

6A 4 w 2 p 2 + A 4 nn A + A 4 mp 4 = 6A 4 /j, 2 u 2 + A 4 n// 4 + A A mv A + 0(rT 1 ^ 2 ). 



In summary, the value of the integrand for (er, t, /i) G Q PI Ai is exp (L 2 + 0(n 1 ^ 2 )) , 
where 

m— 1 n— 1 

L 2 = -Amnp 2 - Anp 2 - Amv 2 - y^ y ( a j* ~ Ax 'j) a ] ~ ^( a *k - A y'k) T k + §A A \i 2 v 2 

j=i k=i 

+ A 4 n/i 4 + A i mv A — iA 3 nfi 3 — iA 3 mv 3 — 3iA 3 cn/j, 1 n 2 — 3iA 3 dmv l v 2 

TO— 1 n— 1 

- 3iA 3 nnfM 2 - 3iA 3 m/iz/ 2 - i /3j»<7- - i ^ P** T l 

3=1 k=l 
rn—l 71—1 m— 1 n— 1 

j,/=i fc,&'=i i =1 fc=1 



with 



n 



2,4m 

V' = ^4^(( X + rfn + d2n2 ) a *k + dna„ k >) = 0(m~ 1/2+£ ), 
Uj k = ^ 3 (naj k + (1 + dn)a^ + cna^) — 3/5 jfc = 0(m~ 1 ^ 2+2e + rT l ^ 2+2e ), 

Vj k = 2^4 3 ( ma i k + + cm ) a *k + dmaj^) — 3(3j k = 0{rrr l ^ 2+2e + rT l ^ 2+2e ). 
Note that the O(-) estimates in the last four lines are uniform over k, k' . 

4.6 Estimating the main part of the integral 

Define E 2 = exp(L 2 ). We have shown that the value of the integrand in Q n M. is 
E l = E 2 (l + 0(n -1 ^ 2 )). Denote the complement of the region M. by A4 C . We can 
approximate our integral as follows: 

E x = f E 2 + 5(n- 1/2 ) I \E 2 \ 
QnM J QnM J QnM 

[ E 2 + 0(n- 1/2 ) [ \E 2 \ 
JQnM Jq 

[ E 2 + 0(1) [ \E 2 \+d(n- 1/2 ) [ \E 2 \. (32) 
JO. Jo.nM c Jo. 



29 



It suffices to estimate the value of each integral in (1321) . This can be done using the same 
calculation as in Section 4.3 of jl], using d jfe = ctj k — Ahj k in place of the variable ctj k used 
in that paper. A potential problem with this analogy is that the variable aj k used in jl] 
has the property aj k = 0{n~ 1 ^ 2 ) 1 whereas it is not true that &j k = 0{n~ 1 ^ 2 ). However, 
a careful look at Section 4.3 of [I] confirms that only the properties = a,* — Ax'j = 
0(n 1//2 ), a^j, = — Ay k = O^ 1 ^ 2 ), and the bounds on g>./, h kk ',Uj k ,Vj k , are required. 

The result is that 



\AmnJ \AnJ \AmJ 



9^ 3A 4 /m n\/3A 4 15 A% 
X ' V|) * 8^ 3 2A 2 Vn m)\AA 2 1QA 3 

n—l 



k=l 

1 ~ \ 



(33) 



where 6 is specified in Theorem 12. 11 

Using ffl3|) and the conditions of Theorem 12.11 we calculate that 



(^ + ^)-^X 2 X + d(n 1 
\ n m J 1 

m—l 

^(^) 2 = i(l-2A) 2 J R 2 + 0(n 3 / 2 ), 



j'=i 

n-l 



^(d, fc ) 2 = i(l-2A) 2 C 2 + 0(n 3 / 2 ). 



fc=i 



Substituting these values into (133|) together with the actual values of A, A 3 ,A 4 , we con- 
clude that 



I" / 7T X 1 / 2 / 7T Y™- 1 )/ 2 / 7T \( 

2 VAmri/ \An) \AmJ 



7T \. 


'" 2 ( 


AmnJ 








x exp| 


a 




+ : 



I - 2 A /m + n \ 1 - 4A /i? > ("'. 



4 + (34) 



24A \n mJ 1QA 2 V 
R 2 + C 2 \ 2 X f l In . 
AAmn A A \m n/ 

By the same argument as in [4j, the other two terms in (132]) have value 0{rT b ) Jq E 2 . 
Multiplying fl34|) by the Jacobians of the transformations T 2 and T 3 , we have proved the 
following. 
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Lemma 4.3. The region S' given by (1301) contains S and 



1 \-2A(m n\ 1-4A/R 2 C 2 



Y <>XP 1 ~2 " "W I n + m) + "W + m 2 

R 2 + C 2 X 2 X f l In 
H r-; h —r-r\ H - + Oin , 



AAmn A A \m n 



4.7 Bounding the remainder of the integral 

In the previous subsection, we estimated the value of the integral I n >(s,t, H), which is 
the same as I(s,t,H) except that it is restricted to a certain region Tl D K. In this 
subsection, we extend this to an estimate of J(s, t, H) by showing that the remainder of 
the region of integration contributes negligibly. 

For 1 < j < m, 1 < k < n, let Aj k = A + a, fc = |A jfc (l — Xj k ) (recall (1231 ). and define 
A min = min jfc Aj k = A + 0{n^ 1 ^ 2 ). We begin with two technical lemmas whose proofs are 
omitted, and a well-known bound of Hoeffding. 

Lemma 4.4. 

if(0»i= n f jk (9 3 +<p k ), 

where 

fjk(z) = \A -4A| fe (l -cos z). 

Moreover, for all real z, 

< fjk(z) < exp(-A jk z 2 + ±A jk z 4 ). □ 
Lemma 4.5. For all c > 0, 

"8tt/75 
t/7 



/8tt/75 
exp(c(— x 2 + ?£ 4 )) dx < yjir/c exp(3/c). □ 
-8tt/75 



Lemma 4.6 ([10]). Let X 1 , . . . ,X N be independent random variables such that EJQ = 
and |Xj| < M for all i. Then, for any t > 0, 



N , 2 \ 

2=1 X ' 
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Lemma 4.7. Let F(0,<f>) be the integrand of I(s,t, H) as defined in (j4j). Then, under 
the conditions of Theorem \2.1[ 

[ \F(0,(j>)\d0d(j) = O(n- 1 ) [ F(9,4>)d9d<f>, 
Jri Jri 

where 1Z C denotes the complement oflZ. 

Proof. Our approach will be to bound J \F(0, cf>)\ over a variety of regions whose union 
covers 1Z C . To make the comparison of these bounds with J n > F(0, (p) easier, we note that 

/ F(0,ci>)d0d(i) = exp(A- 1 O(m £ + n £ ))l = exp(0(m 3£ + n 3£ ))l u (35) 
Jri 



where 



7T X 1 / 2 -i-r/ 7T X 1 / 2 -1-1-/ 7T X 1 / 2 

3=1 J k=l 
m/2 / -jr \n/2 



1 VAn/ Urn/ 



To see this, expand 



A,-, = An + a„ = An exp ( — — + •••), 

3 3 \An 2A 2 n J 

and similarly for A, k , and compare the result to Lemma 14.31 using the assumptions of 
Theorem 12.11 It may help to recall the calculation following fl33|) . 

Take k = 7r/300 and define w , w 1 , . . . , w 299 by w e = 2£k. For any £, let S^i) be 
the set of (0,4>) such that 9j G [wt — k,Wi + k] for at least nm/ir values of j and 
4> k £ [—w £ — 2k, — uy + 2k] for at least n £ values of k. For {0, (f>) G Si(£), 9j + <fi k [— k, k] 
for at least nim — 0(m £ ))n £ fir pairs (j,k) with hj k = so, by Lemma [4.4[ \F(0,<p)\ < 
exp(— CiA min mn £ ) for some c x > which is independent of £. 

Next define S 2 {£) to be the set of (0, <p) such that 9j G [w t — k, w^ + k] for at least nm/it 
values of j, (f) k G [—w^—2k, —w e +2n] for at least n—n £ values of k and 9a ^ 3k, w £ +3k] 
for at least m e values of j. By the same argument with the roles of and 4> reversed, 
\F{0, 4>)\ < exp(— c 2 A min m £ n) for some c 2 > independent of £ when (0, 0) G S 2 {£). 

Now define T^^) to be the set of pairs (0, 4>) such that 9j G [w f — 3k, uy + 3k] for at 
least m — m £ values of j, and <p k G [— — 3k, — + 3k] for at least n — n £ values of k. 
By the pigeonhole principle, for any there is some £ such that [w£ — k, w e + k] contains 
at least nm/n values of 9a. Therefore, 



299 299 

c 



(u 



TlM Q\J{S 1 {£)V}S 2 {£)). 
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Since the total volume of (U£^iW) C i s at most (2m) m+n , we find that for some c 3 > 0, 

/ \F(0,4>)\dOd4> 

< (27r) m+n (exp(-c 3 A min mn e ) + exp(-c 3 A mhl m £ n)) 

< e~ n h. (36) 

We are left with (0, <f>) G U^i(^)- ^ we subtract from each 9j and add -u^ to each 
<p k the integrand F{6, <j>) is unchanged, so we can assume for convenience that £ = and 
that (0,4>) E TZi — TZi(0). The bounds we obtain on parts of the integral we seek to 
reject will be at least 1/300 of the total and thus be of the right order of magnitude. We 
will not mention this point again. 

For a given 6, partition {1, 2, ... , m} into sets J = Jq{0), J\ = J\{0) and J 2 = J 2 (0), 
containing the indices j such that \9j\ < 3k, 3k < \9j\ < 15k and \9j\ > 15k, respectively. 
Similarly partition {1,2, ... ,n} into K = K (<j>), K 1 = ifi(</>) and K 2 = K 2 ((f>). The 
value of \F{0, 4>)\ can now be bounded using 

f jk (9j + 4> k ) 

' exp(-A min (^ + <\> k f + ^Ani„(^ + <Pkf) if U, k) G (J U Jj) x (tf U iQ, 

< < v /l-4A nin (l-cos(12/ t )) < e'^n/ 64 if (j, k) G (J x if 2 ) U ( J 2 x if ), 

1 otherwise. 

Let l2{ m '2i n 2) be the contribution to J n \F{6,<p)\ of those (0,<f>) with |J 2 | = m 2 and 
| if 2 1 = ^2- Recall that | J | > m — m e and |if | > n — n £ . We have 

h{m 2 ,n 2 )< ( m )( n )(2nr^ 

\m 2 J \n 2 J ( 37 ) 

x exp(--^A min (n - 0(n £ ))m 2 - -^ min (m - 0(m £ ))n 2 )l 2 (m 2 ,n 2 ), 

where 

/15k i>15k t j . 

• • • / expf -A min (^ + fc ) 2 + ^A min (^ + fc ) 4 J dO'd^', 
-15k •/ — 15k ^ ., = "" ., rr 

jkeH jkeH 

and the primes denote restriction to j 6 J U J x and G K U K 1 , in the case of the 
summations in addition to the restriction given by the summation limits. Write m = 
m — m 2 and n = n — n 2 and define 6' = (m')^ 1 Y^'j 6j, Qj = Qj — f° r j £ Jo U J 1; 
0' = (n') _1 El 0fc, k = <Pk- 4> for k <E K U K 1} v' = 4>' + 6' and $ = V - </>'. Change 
variables from (0', (/>') to | j G J 3 } U {4> k \ k G if 3 } U {V, ^Z}, where J 3 is some subset 
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of m — 1 elements of J U J x and K 3 is some subset of n — 1 elements of K U 1^. From 
Subsection 14. 31 we know that the Jacobian of this transformation is m'n /2. The integrand 
of ^2 can now be bounded using 

+ ^kf = (n - 0(n £ )) J^9] + (m' - 0(m e )) J^J + (mV - 0{X))v' 2 

jkeH j k 

and 



14 
V . 



+ <PkY < 27n Oj + 27m' ^ + 27m W 

The latter follows from the inequality (x + y + z) 4 < 27(x 4 + y 4 + z 4 ) valid for all x, y, z. 
Therefore, 



0(1) r MR rovn paok / 
l' 2 (m 2 , n 2 ) < ]—f / / ' • • / ex P ( Anin(w' - 0(n £ )) V 0(6/,.) 

m n J -30k J -30k J -30k V 

< 

+ ^min(™' - 0(m £ )) ^(0 fc ) 

+ AutJm'ri - 0{X))g{v'Yj d6 je j 3 d(f> keK3 du, 

where g(z) = —z 2 + \z 4 . Since g(z) < for \z\ < 30k, and we only need an upper bound, 
we can restrict the summations in the integrand to j G J 3 and k G K 3 . The integral now 
separates into m + n — 1 one-dimensional integrals and Lemma 14.51 (by monotonicity) 
gives that 



I 2 (m 2 ,n 2 ) = 0(1)- 



n (m +n )/2 



^(m+n -l)/ 2( y _ _ Q^m' /2-1 

x exp(0(m'/(A mhl n) +n'/(A min m))). 
Applying (I35I) and (|37|) . we find that 



m n 



£ E ^K, n 2 ) = 0(e- C4Am + e~ c * An ) I, (38) 



m 2 =0 ri 2 =0 
m 2 +n 2 >l 



for some c 4 > 0. 



We have now bounded contributions to the integral of \F(0, cf))\ from everywhere 
outside the union of 300 equivalent translates of X — 71, where 

X={(0,4>)\ \0j\, \(j) k \ < 15/« for 1 < j < m, 1 < k < n }. 
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By Lemma [4.41 we have for (6, <fi) £ (j m + n ( w hich includes X) that 

m n 

\F(0, cf>) I < exp (- A ik0i + h + »? + J2 E E + ^ + z/ ) 4 ) ' 

jfcei? i =1 fc=1 

where 9a = 9 3 — 0, 4> k = 4> k — and v = + <f>. As before, the integrand is independent 
of if) = — (f) and our notation will tend to ignore if) for that reason; for our bounds it 
will suffice to remember that if) has a bounded range. 

We proceed by exactly diagonalizing the (m+n + l)-dimensional quadratic form. Since 

EJLi Qj = ELi 0* = °' we have 

m n 

jk&H j =1 fc=l 

m n 

+ 2 ^2 z2( a jk - A jkh jk )9j(j)k 

3=1 k=l 

m n 

+ 2z/ ~ A 3»\h)9j + 2 " E^* fc ~ A »k\H)(f>k- 

3=1 k=l 

This is almost diagonal, because a jk = 0(n~ 1 ^ 2 ), Aj,\ H = 0(1), A. k \ H = 0(1). The 
coefficients —2Aj k hj k can be larger but only in the 0(n) places where hj k = 1. We 
can make the quadratic form exactly diagonal using the slight additional transformation 
(J + U~ 1 Y)~ 1 ^ 2 described by Lemma 14.14 where U is a diagonal matrix with diagonal 
entries Aj t ,fj, A 9k ^ and A t9 ,g. The matrix Y has zero diagonal and other entries of 
magnitude 0(n~ 1 ^ 2 ) apart from the row and column indexed by v, which have entries 
of magnitude 0(n 1 ^ 2 ), and the 0(n) just-mentioned entries of order 0(1). By the same 
argument as used in Subsection 14. 5[ all eigenvalues of U^Y have magnitude 0(n~ 1 ^ 2 ), 
so the transformation is well-defined. The new variables {$j}, {0k} an d v are related to 
the old by 

(§ 1 ,...,e m ,(j) 1 ,..., <j> n , v) T = (i + u^ 1 Yy 1/2 (^ 1 , ...,'& m ,<p 1 ,..., <p n , v) T . 

We will keep the variable if) as a variable of integration but, as noted before, our notation 
will generally ignore it. 

More explicitly, for some d 1 , . . . , d m , d[, . . . , d' n = 0(n~ 3 ^ 2 ), we have uniformly over 



35 



j = 1, . . . , m, k = 1, . . . , n that 

m n 

h = h + E 0(n~ 2 % + d (n- 3/2 + n x h 3k )(p k + 6V 1/2 )z>, 

q=l k=l 
m n 

k = 0k + Yl 0(n- 3/2 + n-% k )^ + O(n- 2 )0 q + 0(n" 1/2 )z>, (39) 

3=1 5=1 
m n 

v = i> + ^ dj^j + ^ dfc^fc + 0{rT l )v. 

3=1 k=l 

Note that the expressions O(-) in (I39j) represent values that depend on m,n, s,t but not 
on {^},{<£j,z>. 

The region of integration X is (m+n) -dimensional. In place of the variables (6, <fi) 
we can use (0, 0, v, ip) by applying the identities 8 m = — Y^^i @j an d 4>n — ~ Ylk=i 



(Recall that 9 and (ft don't include 6 m and <fi n .) The additional transformation ( 1391) 
maps the two just-mentioned identities into identities that define d m and (p n in terms of 
<£, z>), where t? = (i? x , . . . , $ m _i) and (p = (0 l7 . . . , </? n _i). These have the form 

m— 1 n— 1 

(40) 

m— 1 n— 1 

= £ 0(n- 1/2 )^. + d (n- l ))0 k + 0{n'l 2 )v. 

j=i k=i 

Therefore, we can now integrate over (■d,ip,i',ip). The Jacobian of the transformation 
from (6, (f)) to (6, </>, v, ip) is mn/2. 

Next consider the transformation T 4 (i9, (p, v) = (6, <f), v) defined by (139]) . The matrix 
of partial derivatives can be obtained by substituting (HO"]) into (139]) . Without loss of 
generality, we can suppose that x m ,y n = 0(1). Recall that the Frobenius norm of a 
matrix is the square root of the sum of squares of absolute values of the entries. After 
multiplying by n 1 ^ 2 the row indexed by v and dividing by n 1//2 the column indexed by 
v (these two operations together not changing the determinant), the Frobenius norm of 
the matrix is 0{n~ 1 ^ 2 ). Since the Frobenius norm bounds the eigenvalues, we can apply 
Lemma [4.21 to find that the Jacobian of this transformation is 1 + 0{n^ 1 ^ 2 ). 

The transformation T 4 changes the region of integration only by a factor 1 + 0{n^ 1 ^ 2 ) 
in each direction, since the inverse of (159]) has exactly the same form except that the 
constants {dj}, {d' k }, while still of magnitude 0(n~ 3//2 ), may be different. Therefore, the 
image of region X lies inside the region 

y={(#, ip, z>) | 1^1, \<p k \ < 31k (1 < j < to, 1 < k < n), \u\ < 31k }. 
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We next bound the value of the integrand in y. By repeated application of the 
i 2,1 

2 X ' 2< 



inequality xy < \x 2 + \y 2 , we find that 



m n m n 

T2 E E + ^ + ^) 4 < S ( E A ^ + E A -^t + a..v") , 

j=l k=l j=l k=l 

where we have chosen || as a convenient value greater than | (to cover the small variations 
in the coefficients) and less than | (to allow us to use Lemma S3]). Now define h(z) = 
-z 2 + f§/. Then, for (0, <p, v) G y, 



\F(6, 4>) | < exp ( A i*\S Kh) + E A 'k\v h (^) + A„\ B h(u 

0=1 k=l 
m—l n— 1 

< exp ( £ /»&) + A »k\R H0k) + A„ ]R h(u)) (41) 

j=l fc=l 

m—l n— 1 

= exp(A.. |/? /i(z>)) J] exp(A i . |# /i( 1 ? i )) JJexp^^/i^)), (42) 
i=i fc=i 

where the second line holds because h(z) < for \z\ < 31 k. 
Define 

W = {(i>,£,v)ey\ |^-|<|n- 1/2+E (l<j<m-l), 

|&| < \m- 1/2+e (l<k< n-1), 
\u\ < |(mn)- 1/2+2£ }, 

w = y- Wo, 

m—l n— 1 

w 2 = { (£, v) e y I |£ dft + ]T < n~ 5/4 }. 

i=i fc=i 

Also define similar regions Wo, W(, W 2 by omitting the variables t? l5 instead of i? TO , (p n 
starting at (l4"Ti) . (Note that without loss of generality we can also assume that x 1; y x = 
0(1).) Using fl3T?|) . we see that T 4 , and the corresponding transformation that omits $1 
and (fix, map 1Z to a superset of W fl W 2 fl W PI W 2 . Therefore, ^ — 1Z is mapped to a 
subset of Wi U (W - W 2 ) U W( U (W - W 2 ) and it will suffice to find a tight bound on 
the integral in each of the four latter regions. 

Denoting the right side of (j42]) by F G yd, (p, z>), Lemma I4T51 gives 

/ F ('d,<p,v)d'3d<pdv = exp(0(m £ + n £ ))l . (43) 
Jy 
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Also note that 

f 31k 



exp(c h(z)) = 0(1) exp(c /i(z )) (44) 

for cq, z > and z = o(l), since h(z) < h(z ) for z < 2 < 31k. By applying ( 1441) to 
each of the factors of (j4"2l) in turn, 

/ F (tf, 0, z>) dbdtpdu = 0{f^ An?E + e- C6A " 2E )/o (45) 

for some c 6 > and so, by (1431) and (1431) . 

F (tf, <£, z>) dftd(pdi> = exp(0(m e + n £ ))/ . 



Applying Lemma H3J twice, once to the variables d^i, . . . , d m ~ii?m-i> d' 1 (f> l , . . . , d' n _ l (p n _ 1 
and once to their negatives, using M = 0(n~ 2 ), N = m + n — 2 and t = n _5//4 , we find 
that 

/ F (tf,(p,v)d$d<pdv = 0(e"™ 1/4 ) / F Q (0,<p,v)d0d$dv 

= O(e- nl/5 )l . (46) 

Finally, parallel computations give the same bounds on the integrals over W[ and 
W - W 2 . 

We have now bounded f\F(0, <f>)\ in regions that together cover the complement of TZ. 
Collecting these bounds from (133]) . (138|) . (j45|) . (j46|) . and the above-mentioned analogues 
of (USD and (1431). we conclude that 



/ \F(0, 0)| d6W0 = 0(e- c ' Am2e + e- C7An2£ )l 
for some c 7 > 0, which Lemma [4.71 by ([35]) . □ 
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